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Abstract 

Recent observations of the J /ijj spectrum produced in e+e" collisions at the T(4S') resonance are 
in conflict with fixed-order calculations using the Non-Relativsitic QCD effective theory (NRQCD). 
One problem is an enhancement in the cross section when the J/ip has maximal energy, due to 
large perturbative corrections (Sudakov logarithms). In a recent paper, the Sudakov logarithms in 
the color-octet contribution were summed by combining NRQCD with the Soft-Collinear Effective 
Theory. However to be consistent, the color-singlet contributions must also be summed in the 
endpoint region which was not done in that paper. In this paper, we sum the leading and next- 
to-leading logarithms in the color-singlet contribution to the J/ij: production cross section. We 
find that the color-singlet cross section is suppressed near endpoint compared to the fixed order 
NRQCD prediction. 
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I. INTRODUCTION 



Bound states of heavy quarks and antiquarks have been of great interest since the dis- 
covery of the J/V' [H- In particular the production of quarkonium is an interesting probe of 
both perturbative and nonperturbative aspects of QCD dynamics. Production requires the 
creation of a heavy QQ pair with energy greater than 2mQ, a scale at which the strong cou- 
pling constant is small enough that perturbation theory can be used. However, hadronization 
probes much smaller mass scales of order mgf ^, where v is the typical velocity of the quarks 
in the quarkonium. For J/V^, mgf^ is numerically of order Aqcd so the production process 
is sensitive to nonperturbative physics as well. 

Many phenomenological problems can be understood well enough by using the Non- 
Relativistic Quantum Chromodynamics (NRQCD) 0, [sj. NRQCD provides a generalized 
factorization theorem that includes nonperturbative corrections to the color-singlet model. 
All infrared divergences can be factored into nonperturbative matrix elements, so that in- 
frared safe calculations of inclusive decay rates are possible 0]. However, there are some 
predictions of NRQCD are in conflict with the data, in particular the predicted polarization 
of J /if) at the Fermilab Tevatron [5, 6] and more recently the production rate of JhJj asso- 
ciated with extra c and c quarks (both inclusive and exclusive) at the B factories 0, 11] ■ In 
particular. Belle reports a large cross section for J /ip produced along with open charm 0], 



a(e+e- ^ J/^X) 



The predicted ratio from leading order color-singlet production mechanisms alone is about 



0.2 [9|, [lOj and a large color-octet contribution makes this ratio even smaller. In addition to 
the inclusive measurements. Belle reports a cross section for exclusive double charmonium 
production which exceeds previous theoretical estimates. Recent attempts to address the 



latter problem can be found in Ref. [11 



The inclusive J /if) production at the B factories is another potential conflict between 
experimental observations and theoretical predictions using NRQCD [li,[3. Leading order 



NRQCD calculations predict that for most of the range of allowed energies prompt J/ip 
production should be dominated by color-singlet production mechanisms, while color-octet 
contributions dominate when the J/ip energy is nearly maximal. Furthermore, as pointed 
out in Ref. [13], color-octet processes predict a dramatically different angular distribution 
for the J/ip. Writing the differential cross section as 

'^'^ 5(p^)[l + A(p^)cos2^], (1) 



dp^ d cos 9 

where is the J /ip momentum and 9 is the angle of the J /if) with respect to the axis defined 
by the e^e~ beams, one finds the color-singlet mechanism gives A{p^) except for large 
p^, where A{p^) becomes large and negative. On the other hand, color-octet production 
predicts A{p^) 1. The significant enhancement of the cross section accompanied by the 
change in angular distribution were proposed as a distinctive signal of color-octet mechanisms 
in Ref. [13] • It was expected that these effects would be confined to J/ip whose momentum 
is within a few hundred MeV of the maximum allowed. 

However, experimental results do not agree with these expectations. The cross section 
data as a function of momentum does not exhibit any enhancement in the bins closest to 
the endpoint. On the other hand, the total cross section measured by the two experiments 
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exceeds predictions based on the color-singlet model alone. The total prompt J/ip cross 
section, which includes feeddown from ip' and Xc states but not from B decays, is measured 
to be atot = 2.52 ± 0.21 ± 0.21 pb by BaBar, while Belle measures atot = 1.47 ± 0.10 ± 0.13 
pb. Estimates of the color- singlet contribution range from 0.4 - 0.9 pb 0, M, H M 



Furthermore, A{p^) is measured to be consistent with 1 (with large errors) for > 2.6 GeV 
(Belle) and > 3.5 GeV (BaBar). 

The NRQCD factorization formalism shows that the differential J/ip cross section can be 
written as 

daie+e" ^ J/^ + X) = Y^ da{e^e- cc[n] + X) {Oi'^) , (2) 

n 

where da is the inclusive cross section for producing a cc pair in a color and angular momen- 
tum state labeled by [n] = In this notation, the spectroscopic notation for angular 
momentum quantum numbers is standard and z = 1(8) for color-singlet (octet) production 
matrix elements. The short-distance coefficients are calculable in a perturbation series in 
as- The long-distance matrix elements {O'n^) are vacuum matrix elements of four-fermion 
operators in NRQCD 0. These matrix elements scale as some power of the relative velocity 
f <C 1 of the c and c quarks as given by the NRQCD power-counting rules. 

At lowest order in v the only term in Eq. ([2]) is the color-singlet contribution, [n] = ^S^i\ 
which scales as v^. The coefficient for this contribution starts at O(a^) [l3]. There are two 
different contributions to the leading-order color-singlet, depending on what else is produced 
along with the J/ip: e^e" — * J/ip + g + g and e^e~ — J/ip + c + c. Away from the kinematic 
endpoint E^ax = + M^) / {2^/s) , where s is the center-of-mass energy squared, color-octet 
contributions also start at 0(0;^). Since the color-octet contributions are suppressed by 

~ 0.1 relative to the leading color-singlet contributions, they are negligible throughout 
most of the allowed phase-space at leading order in perturbation theory. 

The theoretical situation becomes more interesting, however, near the endpoint. The 
lowest-order, color-singlet term approaches a constant^ 

hm , , , = — —r-^-^(l + r)\ cos^6' . (3) 

z^idzdcosO 27s2mc ^ ^\^l-r J 

where r = Aml/s, and z = Ecc/E^"'^ with E^"-^ = ^/s{l + r)/2, while the lowest-order, 
color-octet piece is singular (proportional to a delta function). Physically, when the J/ijj 
emerges with close to the maximal energy, it is recoiling against an energetic gluon jet with 
energy of order My but invariant mass of order Mx Aqcd / . The degrees of freedom 
needed to describe this inclusive jet have been integrated out of NRQCD, and thus cannot 
be described by the effective field theory. The effective theory which correctly describes 
this kinematic regime is a combination of NRQCD for the heavy degrees of freedom, and 
the Soft-CoUinear Effective Theory (SCET) [H, [3, 20, 21| for the light energetic degrees 
of freedom. Furthermore, the renormalization group equations of SCET will sum the large 
kinematic preturbative corrections which appear near corners of phase space. 

In a previous paper [22| the combination of NRQCD and SCET was used to sum the large 
kinematic logarithms (Sudakov logarithms) which appear in the color-octet contribution 



^ The J/tpcc contribution goes to zero before the kinematic endpoint, due to the non-zero mass of the charm 
quarks. It therefore does not contribute to the endpoint contribution given in Eq. ([3]). 
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near the endpoint. For the color-octet contribution, there are also large non-perturbative 
contributions at the endpoint [23^ which must also be summed into a non-perturbative shape 
function. Since the shape function is unknown, in Ref. 22|] the shape function was modeled. 
Since it is universal, it is possible that it could be extracted from another process (such as 
J/ip photoproduction jlSl)- With the summation of the perturbative corrections and the 
simple model chosen, a good fit to the data was obtained. 

However, to be consistent, the color-singlet contribution should also be summed in the 
endpoint region. This is the goal of the present paper. The kinematic logarithms in the 
J /ip + c + c color-singlet contribution are small, since the mass of the charm quark acts as 
a cutoff. However, we would expect that the summed J/ip + g + g color-singlet rate would 
be suppressed relative to the unsuppressed rate. This would help alleviate the discrepancy 
with the open charm data. However, we would not expect a very large suppression except 
right near the endpoint, and thus do not expect that this will be a solution to the J/ip + 
open charm question. This will be confirmed in our analysis in this paper. The remainder of 
the paper is organized as follows. In Sec II a factorization theorem for J/ ip production near 
the endpoint is developed. Then in Sec HI the Sudakov logarithms are summed, including 
mixing with the J/ip + q + q final state. In Sec IV the phenomenology of the J/ip production 
is investigated, and finally we conclude in Sec V. A similar treatment of nonperturbative 
and perturbative endpoint corrections to the color-singlet and color-octet contributions in 
the inclusive decay T — X + 7 can be found in Refs. [251 . 26| . and we will rely on some of 



the results from these papers. Similar results have been previously reported in Ref. [2 



II. FACTORIZATION 



In this section, we will derive a factorization theorem for e+e — > J/ip + X near the 
kinematic endpoint, where the rate can be factored into a hard coefficient, a coUinear jet 



function and a ultrasoft shape function. The derivation is quite similar to Refs. [2l|, |22|, |2J, 



25 



26l . |28| |. We begin by briefly reviewing the kinematics of the process in the e~^e~ center 
of mass (COM) frame [221]. In the COM frame, the virtual photon has momentum = 
V^/2(n^ + n'^) with the lightlike vectors defined as = (1,0,0, 1) and n^^ = (1,0,0, — 1). 
The J/ip is moving in the z-direction with four- velocity 



2\x^/s 



M, 



(4) 



-'cc 



Here is the J/ip mass and x = {E^ + p^)/^/s. The cc pair has momentum 
Mv^ = Mv^ + h-^yf ^ where M = 2mc and i"^ is the residual momentum of the cc 
pair inside the J/ip. In the J/ip rest frame, £^ has components of O(Aqcd), which get 
boosted in the COM frame to scaling as n ■ £ ~ M^Aqct)/{x^/s), n ■ £ ~ Xy/sAqcn/M^ 
and i± ~ Aqcd- The momentum of the gluon jets is 



Px 



1 - - + (1 

X 



x)n^ 



(5) 



where x = xM/M^. In the end point region the NRQCD factorization formula breaks down 
because NRQCD does not include appropriate coUinear modes. When 1 — a; ~ Aqcd/^, the 
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jet is no longer highly virtual. Since m\/ E\ ~ Kqqd/M -C 1, the gluon jet is composed of 
energetic particles with small invariant mass which must be included explicitly in the effective 
theory. Hence, a new factorization theorem is needed to handle the end point, which can 
be derived using a combination of NRQCD for the heavy quark degrees of freedom and 



SCET[18|, Il9|, l20|, |21j which includes the collinear physics. 

SCET has collinear degrees of freedom whose momentum scales as n-p ~ Q, n-p ~ A^Q, 
and p"*" ~ AQ, soft degrees of freedom whose momentum scales as A and ultrasoft (usoft) 
degrees of freedom whose momentum scales as A^. Heavy quark fields in SCET are the 
same as in NRQCD when considering quarkonium. For e^e~ — > J /ip + X, Q is of order ^/s, 
while A ~ \/\ — X ~ a/Aqcd7^- To the order we are working, operators will contain usoft, 
collinear quarks and gluons and heavy quark fields. Soft fields do not enter to the order we 
are interested and are neglected. 

We match QCD onto SCET at the scale Q by evaluating matrix elements in QCD at 
the scale Q and expanding in powers of A. Each order in A is reproduced in the effective 
theory by the product of SCET operators and Wilson coefficients. All the dependence on 
the large scale Q shows up in the Wilson coefficients. We must include all SCET operators 
which can contribute to the process under consideration at each order of A. These operators 
must respect the symmetries of the effective theory. For e^e' ^ J /ip + X, the operators 



must be invariant under both collinear and usoft gauge transformations [21[. Lorentz in- 
variance is realized in the effective theory by additional constraints on the operators, called 
reparametrization invariance (RPI) [29]. 

In the collinear sector of SCET there is a collinear fermion field ^„_p, a collinear gluon 
field y4^g(soft modes are ignored), and a collinear Wilson line 



WJx) 



^ exp i-Qs-^n ■ j 

perms ^ ' 



(6) 



The subscripts on the collinear fields are the lightcone direction n^, and the large compo- 
nents of the lightcone momentum [n ■ q,q±)- The operator projects out the momentum 
label [isl], n ■ V^n,p = 'P^n,p = n ■ pC,n,p- In the usoft sector there is a usoft fermion field 
Qus, a usoft gluon field At^^, and a usoft Wilson line Y. Using the transformation properties 



for each of these fields under collinear and usoft gauge transformations |21|, we can build 
invariant operators. The coUinear-gauge invariant field strength is 

G^^ ^ --W^[zV^„ + g,A^^^^,zV: + g^A^jW, (7) 
9s 

where 

^1^': = ^V + r1 + -m . D, (8) 

and iV^ = id^ + Qs^ts usoft covariant derivative. RPI requires the label operators 

and the usoft covariant derivatives, which scale differently with A, to appear in the linear 
combination appearing in i'Df^. The leading piece of G^^ is order A and can be written as 
n^G7 = i[P,B% where 

B^^ = -W\Vl + gs{A^^J,_)W. (9) 
9s 

The subscript ± on indicates that /i must be a perpendicular direction. 
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+ cross terms 



FIG. 1: Matching the production amplitude for e^e — > cc + gg in QCD and SCET. Collinear 
gluons are represented by a spring with a line through it. 

We next construct the operators necessary to describe color-singlet ^Si production at the 
end point. A cc pair in a color-singlet ^5*1 configuration must be accompanied by a colorless 
jet of quarks and gluons. The leading operator must have two gluon field operators to create 
the collinear gluons in the final state. Thus, we should construct the operator out of two 
B±^ fields in color singlet configuration. Taking gauge-invariance into consideration, the only 
operator is 

0,,,{1,'S,) = xlpA ■ ^'^^Tr{Bir%';\P,P^) Bi] . (10) 

At leading order, the coefficient is determined by requiring the SCET matrix element of 
Eq. ffTOl) to reproduce the lowest order QCD diagrams for e"^e~ — > cc + gg, shown in Fig. [H 
Matching at tree level, we obtain 

(1,35) 327recea, r^/ 1-r \ 

where r = Aml/s and g'^ = g^'^ — {n^n^ + n'^n^)/2. We can also have a jet made up of a 
quark- ant iquark pair. Again, taking gauge-invariance into account, the only operator is 

0,,,{l,'Sr) = xlpA ■ a'%^^,,WnTff''\v,V^) Wl^n,p. (12) 
The leading order Wilson coefficient is zero. However, since this operator occurs at the same 



order in A, it can be generated through mixing. Just as in the case of Ref . [26| , the mixing 
is small, and we will neglect this term for now. 

At leading order in the SCET power counting the cross section in the endpoint can be 
expressed in a factored form to all orders in a.^ 

1 2 r 

2E^^ = ^^L^''H^,J d/+5(/+,/i)^(/+- v^(l-x)), (13) 

where J is the collinear jet function, S is the usoft function and H^i, is the hard coefficient. 
We shall now prove this factorization theorem. Using the optical theorem, the production 
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cross section can be written as 

'^E^^ = ^^L^'^5^(0|jt(0)|J/^ + X)(J/^ + X|J,(0)|0)(2vr)V(g-p^-px) 

„2 



e2 



L^^ / Ae-'^-^5^(0|jt(^)|j/^ + X)(J/^ + X|J^(0)|0) 



where the sum includes integration over the phase space of X. The lepton tensor is 

L^" =ptP2+P\p'2-9'''Pi-P2^ (15) 
where pi^2 are the momenta of the electron and positron, respectively, and 

T,, = -i I rfV-^^^$^(0|Ji(y)|J/M/ + X)(J/M/ + X|jt(0)|0). (16) 

The first step is to match the QCD current in Eq. (fT^ to leading order in A, 

r = Y, e-^(*'''-^"/')-^r"^'^^J,^5(a;) , (17) 

where the effective current is 

= ^t(A . ^)^x-p{Tr[525.p_i??]} , (18) 

and given in Eq. flTT]) . Substituting Eq. flTTl) into Eq. flTB]) and using q^^ — Mv^ + 

^/s(l — x)n^/2 gives 

^M^' = 5Z r^,^,5,^r„;35i.7;7'^^''^^' (^, a;', f, /i) , (i9) 



where 



T^eff — 
-l-aa'/SP'SS' — ^ 



/" rf4^e-^v^(i-^)^-^ 5^(0| J^,^,,,(^')l^/^ + ^)(^/^ + XMoo)\0) . (20) 



Next we decouple the usoft gluons in Tcs using the field redefinition [21 



<, = ^4T^^ ^ W^„ = Fl^frt^ (21) 



where the first identity implies the second. The collinear fields with the superscript (0) do 
not interact with usoft fields to lowest order in A. In the color-singlet contribution all usoft 
Wilson lines Y cancel due to the identity Y'^Y = 1. Furthermore, the J/tp does not contain 
any collinear quanta, so using 

J2 I W + Xu){J/^p + Xu\ = alJ2 \Xu){Xu\a^ = a\a^, (22) 
^|X,)(X,| = 1, (23) 
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where a|^a^ projects onto final states containing a J/ip, we can write 

X {Q\{Ti[Bi5^,p_B^]{y)}mB'l5^p_B{mm- 
We can use spin symmetry to simplify the usoft matrix element, 
Af A5(0|xlpO-Yp(z/) a;a^ V^t^^-^_p(0)|0) = 

^-5'^ K^K'; {Q\xiyMy) V^t^^X-p(O)lO). (25) 

Then we can use the identity ^^-'A^A^' = {v^v^' — g^^'), where is the four- velocity of the 
J/ip, to further simplify the result. 

We can define a coUinear jet function from the coUinear matrix element, 

{0\{Ti[B^'6^,p^Bf'']{y)}{Ti[B^6^P_B^]m\0) = (26) 
The jet function, Ji^{k~^,^), is only a function of one component of the usoft momentum. 



k~^, which follows from the coUinear Lagrangian containing only the n ■ d derivative 211]. We 
can also define a usoft function 

S(l+ )^ f c-'^+v- (0|X-pO-Vp(y")4»^^P^'x-p(0)|0) ^27) 

Combining Eqs. (|T9l [20l [26], [27]), we can get the fatorization theorem: 

T,. = H^,j t//+5(/+,/i)^(/+- v^(l-x)), (28) 

with 

H,. - ^iOfeS.Wv^' - g''){gT'9T + 9f9fK,^^,s'^^,^s ■ (29) 

Plugging Eq. fl28]) back into Eq. f]T4l) . proves the result Eq. f[T3|) . 

Changing variables from to 2; = Ecc/E"^"-^ with -E^"^' = -\/i(l + r)/2 and integrating 
over cos 6', we finally get 

= a„Plr, z] j dl+S{l+, ii)J„{l+ - v^(l - £■)) . (30) 



Here 

25677 a'alel (1 + r)(l + 2r) ^^^^^3 ^ 

^0 = ^1 — li::^ ^ — :^ (Ca^ij) (3i) 



is the differential cross section at the end point predicted by NRQCD and P[r, z] 

is a phase space factor. Note that -P[r, 1] = 1. 
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FIG. 2: Feynman diagram for the leading order jet function. 

To leading order the jet function can be calculated easily. The Feynman diagram for the 
vacuum matrix element is shown in Fig. [2l By evaluating the one loop integral, we get 

1 

ImJ^(/c+) = - J . (32) 

Substituting it into differential cross section in Eq. (1501) and summing over u gives, 

^ = aoP[r, z] J /i)0(/+ - v^(l - x)) . (33) 

The color-singlet usoft function just shifts the endpoint from the partonic to the physical 
hadronic endpoint [2^. To show this, we first note that the usoft function can formally be 
written as 

^,^+ . ^ (0|x Wp'^(m ■ d - i+)a\,a^^l(T''x-p\0) 
Then by integrating over in Eq. (1551) gives 



da 

dz 



(ToP[r, z] 



ippQlin ■ d + ^/s{l - x)]a\,aii,ijjl(T^X-p\^) 



4m,(Of(3^i)) 



(35) 



Finally, writing x in terms of z 



X 



sx + Ml/x 
s + Ml 
1 + r 
1 — r 



1 — r 
1 + r 



and using the result in Ref. [30| we get 

da 
dz 



e(l-z)aoP[r, z] 



(36) 
(37) 

(38) 



Notice that as 2; 1, this coincides with the lowest order NRQCD result in the same limit. 



III. RESUMMING SUDAKOV LOGARITHMS 

One of the main strengths of using an effective field theory is the ability to sum logarithms 
using the renormalization group equations (RGEs). Large logarithms of the ratio of well- 
separated scales arise naturally in perturbation theory, which can cause a breakdown of the 
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perturbative expansion. By matching onto an effective theory, the large scale is removed to 
be replaced by a running scale /i. After matching at the high scale, the operators are run to 
the low scale using the RGEs. This sums all large logarithms into an overall factor, and any 
logarithms that arise in the perturbative expansion of the effective theory are of order one. 

For e^e^ — > J/ip + X, there are logarithms of log(l — z) that appear in the perturbation 
series. Near the endpoint, z ^ 1, these become large, and need to be summed, which the 
RGEs will do for us. For the color-singlet ^5*1 contribution, unlike the color-octet process (22| . 
these endpoint lo gari thms are single, not double, logarithms. A similar situation occurs for 
radiative T decay [25|. Double logarithms occur when there is an overlap of soft and collinear 
logs. For the color-singlet case, the soft logarithms do not occur. This can be seen by the 
fact that the usoft Wilson lines canceled out of the color-singlet matrix element. Physically, 
the long-wavelength gluons do not couple to the tightly bound color-singlet cc. 

We have matched in the previous section onto the SCET color-singlet operator, which 
intergrates out the large scale fin, replacing it with a running scale fi. We now run the 
color-singlet operator from the hard scale to the collinear scale, which sums all logarithms of 
1 — z. To run the color-singlet operator given in Eq. flTU]) . we calculate the counterterm for 
the operator, determine the anomalous dimension, and then use this in the RGEs. Luckily, 



the calculation of the anomalous dimension has already been done in Ref. [25|], and we can 



lift the results from that paper. The result for the resummed, differential cross-section is 



dar 



where 7 is defined as 



7 



2 

Jo 



dz 



-11 



aoP[r, z]eil 



df] 



^1—7] 1] 



V))] 



1l 

3 



(39) 



(40) 



To sum the large logarithms, we use the same hard scale as in Ref. [22 

z^H in the above expression.^ 



^H = {slM){l-r) 



and the collinear scale Hc ~ a/F 

To be completely consistent, we should include the mixing of the gg jet with the qq jet. 
Since the match onto the qq operator begins at a higher order than the gg operator, except 



for very close to 2; = 1 the mixing term is small |26l]. The calculation of the mixing in SCET 



was first done in Ref. j26|], and we just quote the results here. Once we included the mixing 



effect, the resummed differential cross section becomes 



1 dar, 



(To dz 



\P[r,z\e{l-z)Y, 



nodd 



qf(") (n)2 



where r(/i) is defined as 



r(/i) 



tts(/iH) 



.(41) 



(42) 



^ The hard scale fXH that we use is different than the choice of Ref. 27[. However, numericaUy they are 
almost the same, and will not have a large effect on the results. 
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and 



(„) _ + l)(n + 2)(n + 3) 



5/2 



An) 
J 3/2 



9(n + 3/2) 
(n + l)(n + 2) 
n + 3/2 



(43) 
(44) 



We also defined A|^^ and 'j^^as 



An) 



(n) 
7± 



(n) 



2 [^as"* ~'~ ^^'J 



±A] 



(n) 
795 



A 



in) 



A 



with 



A 



(n) 
799 



(n) 
799 



(n) 
799 



(n) 

Iqq 



V(7£^-7j?)^ + 47i?7£^ 



+ 



1 1 



+ (n + 2)(ra + 3) 6 z 
1 + 3n + 4 



-nt 



3(n + l)(n + 2) 
+ 3n + 4 



+ l)(n + 2)(n + 3) 



(n + l)(n + 2) 2 



n+l , 

9 ^ 7 



j=2 



In Fig. [31 we plot the difference of the mixing result, Eq. ( HTil . and the non- mixing 
Eq. fl5^ . normalized to the mixing result. For this plot, we chose the scale /ic = VT^ 
The difference between the two is a fraction of a percent, except extremely close 
endpoint, where our results no longer hold. We can therefore use either the mixing 
non- mixing result, Eq. (139!) or Eq. fHTl) . 



(45) 
(46) 

(47) 
(48) 

(49) 
(50) 

(51) 



result, 

to the 
or the 



IV. PHENOMENOLOGY 

The result from the previous section, Eq. (139!) or Eq. (HTl) . summed up the leading log- 
arithmic corrections which are important near the endpoint. Away from the endpoint, the 
logarithms that we have summed are not important and contributions that we neglected in 
the endpoint become important. We therefore would like to interpolate between the lead- 
ing order color-singlet calculation away from the endpoint and the resummed result in the 
endpoint. To do this, we will define the interpolated differential rate as 



1 da: 



int 



(Jo dz 



1 r/rr'^''' 



dz 



-P[r, z] 



1 dcTj-Q^f^ 

ctq dz 



(52) 
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FIG. 3: The difference between mixing and non-mixing da^. 
calculated at the scale /ic = Vl — z^h- 



xjdz, normalized to the mixing result, 



The term in parentheses vanishes as z ^ 1, leaving only the resummed contribution in that 
region.^ Away from the endpoint the resummed contribution combines with the —P[r, z] to 
give higher order in asifJ^n) corrections. 

For our figures, we will use rric = 1.4 GeV and ^/s = 10.58 GeV. In Fig. HI we compare the 
resummed, interpolated result, Eq. (132]) . to the leading-order e"'"e~ — >• J/ipgg color-singlet 
result [9]. We also show the scale dependence of the interpolated result. The dot-dashed 
curve corresponds to the leading-order color-singlet result. All curves are normalized to 
ctq given in Eq. fl3T]) . The solid curve is the interpolated result, plotted at a scale /ic = 
^(1 — z)^H. The dashed curve is the interpolated result at a scale /ic = 2a/ (1 — z)hh, 
while the dotted curve uses the scale fXc = a/(1 — z)fiH/2. As can be seen, there is not a 
large scale dependence. 

As shown in Fig. HJ the resummed result is smaller than the leading order result. In order 
to better see the effects of the resummation, in Fig. [5l we plot the difference of the leading- 



^ This choice of interpolating between the results is different than the one made in Ref . 27| . Given the fact 
that the function P[r, z] is a phase-space factor, we believe our choice more accurately encompasses the 
deviation due to higher-order QCD corrections. The choice in interpolating factor is the largest difference 
between our result and the result of Ref. 27| . Note that the choice made in Eq. ([5^ switches from the 



leading-order result to the resummed result closer to the endpoint than the choice in Ref. 27 1 
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FIG. 4: The color-singlet differential cross section . The dot-dashed curve is the leading-order 
NRQCD preciction. The solid curve is the interpolated result, Eq. [5^) prediction at calculated 
at the scale fic = y^(1 — z)fiH- The dashed curve is the interpolated result at the scale fic = 
2\J (1 — z)iJ,H, and the dotted curve is the interpolated result using the scale fic = — z)hh/2- 



order, color-singlet result and the interpolated result, normalized to the leading-order result. 
As can be seen, in the endpoint region there corrections become large. However, over most 
of phase space, the corrections are less than 10%. 

The total color-singlet contribution also has the J/ip -\- c -\- c final state, so we need to 
combine the results above with the color-singlet e~^e~ — J/ip + cc contribution [9]. In 
Fig. E] we compare the total leading-order, color-singlet result (dotted line) to the total, 
resummed color-singlet result (solid line) for {l/aQ)da/dp^. Also shown as the dashed line is 
the J/ip + c-\-c contribution. While the resummed result is slightly suppressed compared to 
the leading-order result, qualitatively the plots are the same. Note that this implies that the 
resummation of the color-singlet contribution are not big enough to explain the anomalously 
large contribution to J/ ip associated with extra cc found at the B factories 0, [1] . 

In Fig. [3, we plot the color-singlet prediction for A{p^). The dashed curve is the leading 
order, color-singlet result, and the solid curve is the interpolated result, including the J/'ip + 
c + c contribution. Since the resummation is independent of the angle, both curves drop to 
the same value at the endpoint. 



ml 



^(pD = -rvdr- (53) 



s + m. 
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FIG. 5: The difference of the leading-order NRQCD e^ e~ J/tpgg differential cross section and 
the interpolated result, Eq. \5S\) . normalized to the leading-order result. The interpolated result was 
calculated at the scale fic = y^(1 — z)^h- 



Away from the endpoint, the resummed color-singlet rate is slightly larger than the leading- 
order rate. However, to explain the data, we still need to include the color-octet contribution. 
To make a prediction for the differential cross section, we need to combine the color- 



singlet results discussed in this paper with the resummed color-octet results from Ref. [22 
Given the size of the corrections found in this paper, the results are qualitatively the same 
as those presented in Ref. 13 . 



V. CONCLUSION 

In this paper, we study the color-singlet contribution to J /ip production in e+e~ collision 
near the kinematic end point by using a combination of SCET and NRQCD. The calculation 
consists of matching onto a color-singlet operator in SCET which integrates out the hard 
scale. By decoupling the usoft modes from the collinear modes using a field redefinition, we 
are able to show a factorization theorem for the differential cross section. The differential 
rate can be factorized into a hard piece, a collinear jet function, and an usoft function. As 
pointed out by Ref. |30] the usoft function in this case can be calculated, resulting in just a 
shift from the partonic to the physical endpoint. 

By running the resulting rate from the hard scale to the collinear scale, we sum the 
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FIG. 6: Comparison of the leading-order and resummed total color-singlet results. The dashed 
curve is the NRQCD prediction for e'^e" J/tpcc. The dotted line is the total leading-order, 
color-singlet NRQCD prediction, while the solid curve is the total color-singlet prediction including 
the interpolated e'^e~ — > J/i^gg result. The resummed result was calculated at the scale //c = 
1/(1 - z)fiH. 

logarithms of the ratio of the hard and colhear scales, which correspond to large Sudakov 
logarithms of 1 — 2. Finally, we combine the SCET calculation with the leading order, 
color-singlet NRQCD result to make a prediction for the color-singlet contribution to the 
differential cross section over the entire allowed kinematic range. If we combine the results 
for the color-singlet calculation given in this paper the resummed results for the color-octet 
calculation given in Ref. j22|, we now have a consistent prediction over the entire kinematic 
range for the e~^e~ — > J/ip + X differential cross section. 

To be consistent the resummation of the color-singlet presented here must be included. 
However, except for right near the endpoint the size of the corrections are small. The color- 
octet contributions, as can be seen from Ref. [22|, are necessary to get a reasonable fit to the 
data and are larger than the color-singlet contribution over all of phase space. Therefore, 
while the quantitative picture changes slightly, the qualitative picture is the same with or 
without running as what was presented in Ref. j22|. In particular, we still do not have an 
explanation for the unexpectedly large number of J/ip being produced with extra charm. 
The solution to this puzzle will have to come from another source. 
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FIG. 7: The color-singlet contribution to A{p^). The solid curve is the SCET prediction, with 
fic = — z)iJ,H and the dashed curve is the lowest-order NRQCD prediction. 
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